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Entropy Changes in Rarefaction Waves ] 

Robert F. Dressier 

Frictional flow of a polytropic gas is investigated for centered rarefaction waves. En- 
ergy balance is maintained by assuming that mechanical energy loss due to ("fictional force 1 
reappears as heat. Entropy behavior at the wave front is discussed. The aonhomogeneous 
linear system of partial differential equations with variable coefficients for first-order frac- 
tional effects is derived. Utilizing the geometrical similarity of the Mach lines, the func- 
t ional form of these quantities is ascertained, which permits explicit solutions of t he boundary 
value problem. The expansion procedure is singular at the wave front, and results are not 
applicable there. First-order effects are expressed as polynomials plus terms singular a1 
the wave front. Results are compared with expressions obtained when heat generation 
due to frictional force is neglected. 

1. Introduction 

Equations for one dimensional, unsteady flow of a compressible gas permit explicit solu- 
tion for centered rarefaction waves, such as occur in shock tubes or pneumatic control circuits. 
For some applications, however, it becomes important to include also the complications re- 
sulting from frictional dissipation. 'Plus matter is discussed here in two ways, first with the 
simplifying assumption that the frictional effects consist merely in a retarding force, and 
second, by the more complete consideration that includes the resulting heat generation and 
accompanying change in entropy. This means that the mechanical energy destroyed by the 
retardation force reappears as heat energy, resulting in additional expansion. The simpler 
flow description without heat generation will be referred to as model A, and the second, more 
accurate, one as model B. The equations for steady frictional flow in model B have been 
studied and yield energy theorems and an extended form of Bernoulli's theorem [l], 2 but knowl- 
edge of the unsteady flows appears to be restricted mainly to a few isolated numerical calcu- 
lations. 

A centered rarefaction wave is perhaps the most basic pattern in unsteady flow; the present 
paper derives solutions for first-order dissipative effects in this unsteady wave using the lull 
system of equations of model B, and compares results with corresponding ones for model A 
(which were previously published by the author [2] as a supplementary result in a study of 
shock tubes with varying cross section). 

A recent paper by Ludford and Martin [3] lias discussed anisentropic effects in centered 
simple waves. They consider flows where the specific entropy varies from particle to particle, 
but remains constant for a given particle. In frictional flow, however, the specific entropy 
for a given particle must vary with time, and a centered rarefaction wave cannot remain a 
simple wave. 

Consider a polytropic gas at rest for time t <0 in a duct of uniform cross section, filling 
the duct to the left of a diaphragm wdiere x<^0, and with a vacuum initially at 7>(). Heat- 
exchange between the gas and the duct, heat conduction and radiation, and viscous effects are 
disregarded. When the diaphragm is removed at t = 0, a centered rarefaction wave begins. 
Dimensional quantities are denoted by bars in order to define the unbarred dimensionless 
quantities to be used subsequently. Before flow, the gas state is described by sound speed 
c , mass density p„, and temperature T ; after flow we have velocity u, local sound speed c, 
mass density p, and pressure J). 

For model A, assuming no internal heat generation and no entropy changes, the equation 
of state can be written in forms 

p=A*^\ c 2 = y A*^- l =yRT, (1) 



1 This research was supported hy the Office ot Naval Research, CJSN. 

2 Figures in brackets indicate the literature references at the end of this paper. 
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where 7 is the adiabatic exponent, R the gas constant, T the temperature, and A* is a known 
constant (consistent with the assumption here of no entropy change). In this simpler model, 
the friction produces merely a retarding force per unit mass of magnitude \u 2 where A is the 
known friction factor. There is some question concerning the constancy of this coefficient for 
frictional flows that are highly unsteady; this question is discussed by Schultz-Grunow [4], 
Jenny [5], and others. Here it is considered to be a constant. 

Introducing dimensionless (unbarred) variables defined by x=x/D, t=(c,jD)t J u=u/c , 
c=c[c , p^^l~p 0j \=D\ where D is some length associated with duct cross section, the momen- 
tum and continuity equations for model A are 



u t +uu x + 



7-1 



■ cc T =—\u\u\ 



2 2 

7—1 7—1 



(2) 



for the two unknowns u(x,t) and c(x,t); and the energy equation is neglected. 

In deriving equations for model B one must take account of the mechanical energy lost 
through the action of the retarding force \u\u\. The action of a body force cannot itself di- 
rectly change entropy; however, if a term is included in the energy equation for an amount 
of generated heat equal to the lost mechanical energy, this heat will change the specific en- 
tropy. In this way, the energy equation becomes expressible directly in terms of the friction 
coefficient. The equation of state for model B now includes the entropy dependence, 



p=A(!)t\ A(s) = (y-l)e {s - s « )/c » ] 



~-yRTQi,T) 



(3) 



where c v is specific heat at constant volume, s is the specific entropy (entropy per unit mass) 
and lo is a constant with a value depending unon the entropy level adopted. Letting dimen- 
sionless specific entropy s be defined by s=( r J 1 /~c 2 )s, the momentum, continuity, and energy 
equations for model B become, respectively, 



u t +uu x 



7-1 



cc x =—\u\u\ J r c\ 



CU T 



7-1 



uc x - 



7—1 C 



s t J r us x =\ 



(4) 



for the three unknowns u(x,t), c(Xjt), and s(x,t). We set s=s=0 for the gas when initially at 
rest, and the equality 

A* = A(0) = (y-l)e-°^ (5) 

determines the constant ~s . 

2. Characteristic Equations 

The Mach lines for model B are defined by the system of characteristic equations equivalent 
to the totally hyperbolic system (4). These consist of three families of real curves with direc- 
tions given by 

dx 



dx . 
ft: Tt =u+c 

dx 



(6) 
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Figure I 



which ;uc respectively the backward sound path, the forward sound path, and the particle 
trajectory (sec fi^. 1). The corresponding' differential relations along each family are 

ft: d(u— —c]+cds=—\u\u\(l+[y—l]-)dt 



ft: 



<fe=Xi^-<ft. 



(7) 



Those equations arc used in [6] as the basis for a discussion of a numerical computation on a 
rarefaction wave in a shock tube. 

3. Entropy Behavior at the Wave Front 

From the above equations and the equation of state (3), the behavior of s at the forward 
wave front where p=0 can be inferred, consistent with the assumed model for the flow. For 
notational simplicity and to permit numerical results, we now specialize the gas to be air with 
7 = 7/5; however, the same arguments will apply and analogous solutions can be derived for 
any other admissible value of 7. 

lf~So is chosen as indicated and the relation 7=1 -\-R/c VJ is used, we obtain the dimensionless 



equation of state in the form, c 2 



«7(7-l)s -7-1 



p y \ which for air is 



p = e 7s/ V. 



(8) 



On the wave front curve where p=0, this permits two possibilities: either (a), s remains finite 
and hence c = there, or (b), s becomes infinite and c is not known in advance. If s were finite 
and continuous and c vanished there, then by the third equation in (7), the rate of change of s 
along the wave front would be everywhere infinite asu^O there. This follows from the fact 
that the wave-front locus is also a particle trajectory and therefore a f 3 Mach line. Therefore 
s must actually be infinite on the wave-front curve, and c is not determined there by this argu- 
ment. Because s is entropy per mass, this does not necessarily imply that entropy per volume 
becomes infinite at the wave front for this model. 



4. First-Order Effects 



By setting up asymptotic expansions of form 

v(x y t,\)~v\x,t)+V A (x,t)\+ . 
v(x,t,\)~vXx,t)+V B (x,t)\+ . 
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. . (model A) 
. . (model B) 



(8) 



a perturbation scheme will be denned to give information for small values of A. The quantities 
v(x,t,\) represent each unknown u,c,s, and their various first derivatives. The unperturbed 
centered rarefaction wave is the known simple-wave solution 



u°=5(l+x/t)') 



^0=1^ — 



i(5-x/t) 



s°=0 



(10) 



which when \=0 satisfies system (2) for model A and system (4) for model B. The first-order 
terms for model A are U A , C A , (and S A =0) and have been published in [2]. Now we proceed 
to solve the more complicated problem for U B , C B , and S B and to compare these with the 
previous result's. 

The expansion s^0 J rS B (x,t)\-\- . . . shows at once that the perturbation procedure for 
model B cannot be expected to retain validity in the neighborhood of the forward wave front, 
since s(;x,t,\)=co for any A>0 on the wave front whereas s(x,t,0)=0 there. The results 
therefore will not apply in the neighborhood of the forward wave front where a type of boundary 
layer effect occurs. A separate analysis would be required in that region, possibly similar to 
the application of the Pohlhausen method as used by Whitham [7] for a water wave problem.. 
There would be little justification for it in this problem, however, because the model assumed 
for frictional dissipation would in any case become inaccurate in this region of extremely low 
density. (A singular situation does not occur in the perturbation for model A, In that case 
there is no edge layer effect, and results hold throughout the full wave; but on the other hand, 
the basic model itself is everywhere less accurate than model B.) 

When expansions of type (9) are inserted into system (4), after using (10), the first-order 
terms satisfy the system: 



(5t-x) 2 30(t + x) 30(5t-x± 30 rTB , 30 rm , 25(t + x) 2 ' 

5(<+x) obi fiOB I25(t + xf 



(11) 



The solution of these equations that is sought must satisfy the conditions LT B = C B = S B =0 on 
the straight line t=—x, (t>0). This nonhomogeneous linear system has the same character- 
istic curves as the solution u°, c° plus its particle trajectories. By studying these curves one 
can infer an easy procedure for obtaining the desired functions. These curves in the flow 
wedge are 



J*?: x= m t, — l<m<5 

* a)-©-©"" »<« 



(12) 



where each equation defines a one-parameter family of curves over the range of a parameter as 
indicated. The f° and f§ curves emanate from the back wave line x=-— t at a point (— a } a) 
for any a>0 (see fig. 2). Writing {*} as (xja)=m(t/a), one sees that the coordinates of an 
intersection point of a fixed f? curve with any f§ curve must be proportional to the parameter a. 
Therefore as the parameter a is varied, the quantities t/a and x/a at such intersection jDoints 
remain constant as we progress along a fixed f? ray. Likewise an analogous situation will 
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Figure 2. 



apply for intersections of f§ curves with a fixed f? ray. This geometric similarity of the char- 
acteristics can be employed to find the functional form of the desired solutions (7 B , C B , S B . 
The characteristic relations associated with (12) are 



*?: d(l 



' IJ -5r 1! )+( r> ' w -' r )rfs B =( 



25(.r+0 3 25(x+t) 2 



'3 



\dt---—— it 
'-(x—tyT) ,«>r / A t 

jo. ^r/B,^B W A f/ o., / 25(/+<) 3 25(:r+f) 2 \ 

H- rtv ' +0( V 6* /> ' Vl-8^ 2 (^-5<) + 36* 2 J 

,«. ,qb 125(f+x) 3 



icft 



(13) 



The last equation does not contain unknowns on the right side and could therefore l>c integrated 
along any fjj curve. If this were done, using the third equation in (12) and the geometric 
similarity, it is easily seen that it would lead to a relation of the form S B /a = F :i (t/a) when 
S B = at the origin. Because t/a is constant along a fj ray, S B is proportional to a, and hence 
to / there; the functional form must then be S B (x, t)=L(in)t. Using this result for dS B in the 
second equation of (13), by analogous argument applied to the & curves, it follows that U BJ r 
5C B =F 2 (m)t. Finally, we see that if we integrate the first equation of (13) along any .<*',' i'ay 
(where m== constant), we would have (ZT 3 *— 5( m ) proportional to t if the quantity 5{U B — G B )/3t 
were a constant along the ray. Then it would follow that U B and ( m woidd separately be pro- 
portional to t also. But in that case, 5([/ B — C B )/3t would actually be constant as desired. 
There these considerations indicate the existence of a solution to (13) and (11), vanishing at the 
origin as required, with the functional form j(m)t for each unknown U B , ( 1B , and S B . This 
information now permits immediate solution of the original system (11), as follows: 

First, system (11) is transformed through use of the new independent variables m and /, 
putting F B {x, £) = £7*(m, t), and likewise for the other unknowns. Then setting U*=h(m)t, 
C* = k(m)t, and S* = L(m)t leads immediately to a nonhomogeneous system of three ordinary 
differential equations for h, k, and L. Although this system contains variable coefficients, it 
can nevertheless be solved exactly. The particular solution satisfying the conditions h{— 1) = 
Jc(— 1)=L(— 1)=0 is then finally obtained to give the first-order effects for model B. 






7(l + m) 



1_868 , 90_ 

7 + (5-m) ' (5— m) 



] 



S B 

t 



i2- r (5 ~ m)6 *(5 

"^ 5 L2,177,280 5 [ ° 



-m)+3- 



108 



27 



7(5- m) ' (5-m) 2 
269 



] 



M14) 




8 




1 ' 


1 / 


1 


6 


- 




/TOoo 






MODEL 


B \ / 




4 


- 






" 


2 


_ 








** 










2 

O 
o 

<r -2 
o 
















" 


-4 
g-6 










> 




MODEL 


A -^ \ 




-8 








" 


-10 








\ 


-12 








\ - 


-14 




1 1 , 




\ 



2 3 

m( = x/t) 



FlCURE 4. 



For comparison of these results with the corresponding expressions for the simpler flow 
model, we list below the expressions obtained in [2] for model A: 

U A 25f2^ , 3 K , R 

— =— r -(5— mY— 5(b— 



^| |~| (5-m) 3 -5(5-m) 2 +36(5-m)-84l 
*Y=Y^ [~ (5-m) 3 -10(5-m) 2 +72(5-m)-168+3(m + l) 2 l 



=0. 



(15) 



The expressions in (14) have singularities at m = 5, at the forward wave front for the unperturbed 
flow. This is to be expected in view of the previous remarks concerning the entropy behavior 
at the forward wave front that makes the perturbation procedure singular. On the other 
hand, the expansion for model A does not possess this singular property, and the expressions 
(15) do not become infinite at m=5. A comparison of (14) with (15) must be restricted to a 
region that excludes the forward wave-front zone, Figure 3 gives the numerical evaluation of 
the quantity S B /t. When drawn on a semilogarithmic plot, this curve is almost linear over the 
middle range l<ra<4, and for convenience can be accurately approximated there by the 
expression .63 exp (1.95 m). Figure 4 compares the behavior of U B /t with CJ A /t. In the first 
half of the flow zone, the two models predict almost identical results, but further on the velocity 
becomes higher for model B. In equations (4), as s z >0 for this case, the effect of the entropy 
term in the momentum equation is to diminish the total resistive force, and this by itself would 
increase the velocity. In the continuity equation, however, the extra positive term on the 
right side would have the effect by itself of decreasing the velocity in the back portion of the 
flow zone, because this term is analogous to the term that would be present for a converging 
duct in quasi-one-dimensional flow, as derived in [2]. Evidently the fact that the velocity 
term for model B is slightly below the model A curve in the beginning of the flow zone, is 
due to the extra term in the continuity equation dominating over the entropy term in the 
momentum equation until s z grows sufficiently large to raise the velocity. 
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Figure 5. 

Tho first-order corrections for the local sound speed c are shown in figure 5 for the two 
models. In this case there is no resemblance in (he behavior of the two curves, as the sound 
speed, which is proportional to the square root of the temperature, is more directly sensitive 
than the velocity to the internal heat generation assumed in model B. 

Finally we use our results for sound speed and entropy in the two models to calculate the 
corrections for the dimensionless density p=~p/'p () . The dimensionless equation of state of air 
for model A corresponding to (1) is p = c 5 . After inserting expansions here, the coefficient of X 
in the expansion of type (9) for p is then equal to 5(c°)Y TA . This quantity, after division by /, 
is shown in figure 6. For model B, when corresponding expansions are inserted in (8), the 
coefficient of X in the p expansion is equal to 5(c°) 4 [C B —(7c°S B /25)]. This quantity, after 
division by t, is likewise shown in figure 6, and one observes that the widely different results 
of the two models for sound speed and entropy nevertheless combine to produce density cor- 
rections that are in close agreement. The pressure and temperature behavior for the wave can 
be obtained directly from (3) using the results already presented. 
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